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The pseudoscalar sum rules of the heavy-light quark systems are used for extracting simultaneously the c, b 
quark masses and the decay constants /d (3) ,/b (s) of the D^,B^ mesons. To order a s , one obtains the 
running quark masses: fh c (m c ) = (1.10 ± 0.04) GeV, fh^mi,) = (4.05 ± 0.06) GeV, the perturbative pole 
masses: M c = (1.46 ± 0.04) GeV, M b = (4.69 ± 0.06) GeV, and the decay constants: f D = (205 ± 20) 
MeV, f B = (203 ± 23) MeV and f Ds = (235 ± 24) MeV , f Bs = (236 ± 30) MeV, in the normalization 
where f n = 130.56 MeV. The fitted values of the pole and running masses satisfy quite well their three-loop 
, perturbative relation. The value fo — Ib confirms earlier findings from the sum rule that the 1/ '\/Mp heavy 

qq ' quark symmetry scaling law is affected by large 1/Mp corrections. 
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1 r-| 1 Introduction 

JL , One of the most important parameters of the standard model is the quark masses. However, contrary to 

the leptons, where the physical mass can be identified with the pole of the propagator, the quark masses are 
difficult to define because of confinement. Some attempts have been made in order to define the heavy quark 
pole mass within perturbation theory, where it has been shown to be IR-finite [Q and independent of the 
choice of the regularization and renormalization schemes used Q . More recently, it has been noticed, in the 
limit of a large number of flavours, that the resummation of perturbative series can induce a non-perturbative 
term, which can affect the truncated perturbative result, and can, then, limit the accuracy of the pole mass 
determination (for reviews see e.g. Q|). One may bypass the previous problems, by working, at a given 
order of perturbative QCD, with the running quark masses, which are treated like coupling constants of 
the QCD Lagrangian (see e.g. ||), and where non-perturbative-like effect is expected to be absent. It is 
also known that the decay constants Jd,b of the pseudoscalar D, B mesons play an essential role in the 
neutral D-D and B-B oscillations, and in D, B decays, and can directly be measured in the D~ , B~ — > lv\ 
leptonic decays. In addition, it is important for your knowledge of heavy quark symmetry, how the value of 
these decay constants deviate from the naive 1 / \JMp scaling law, expected to occur when the pseudoscalar 
meson mass is infinitely large. A lot of efforts has been furnished in the literature H for extracting directly 
from the data the running masses of the light and heavy quarks and the heavy quark "perturbative" pole 
masses using the SVZ QCD spectral sum rules (QSSR) (?[ (for a complete review, see e.g.: whilst fp>,B 
come from different forms of the pseudoscalar sum rules Q-fl5| since the pioneering work of |l(J. In this 
note, I shall consider a direct extraction of the running charm and bottom quark masses from pseudoscalar 
two-point function sum rules where the a 2 s correction has been recently obtained in Jljfl and where we shall 
also use the observed values of the meson masses: M D - — 1.869 GeV and M B - — 5.279 GeV. In the same 
time, we shall simultaneously estimate the decay constants /d (s) , /b (s) of the D( s ) and B^ mesons. The 
simultaneous extraction of the quark mass and fp together with the extraction of fp using the running 
quark mass has been initiated in our previous work [ pT| fj"3| to order a s , and will be improved here to order 
a 2 s . All previous works [||-|l^| with the exception of the recent work in |l5| have used, as input, the pole 
mass value for extracting /p, where, as we have mentioned previously, the definition of the pole mass might 
be affected by some non-perturbative contributions Q|. Moreover, the extraction of the quark mass value 
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from the pseudoscalar sum rule itself which is an improvement of our previous works []l4| , |l3f is not done in 
all previous works. 



2 The QCD spectral sum rules 

We shall work with the pseudoscalar two-point correlator: 

Wl 2 ) =ijd 4 x e*** (0\TJ q (x)jl(0)\0), (1) 

built from the heavy-light quark current: Jd(x) — {itlq + md)Q(i"/5)d, and which has the quantum numbers 
of the D and B mesons. rriQ is the heavy quark mass, and we shall neglect the d quark mass here. The 
corresponding Laplace transform sum rules are: 



r°° i d 

CM= dt e- tT - W> 5 (t), and TZ(t) = -— log£(r), 
Jt< T « T 



(2) 



where t< is the hadronic threshold. The latter sum rule, or its slight modification, is useful, as it is equal to 
the resonance mass squared, in the simple duality ansatz parametrization of the spectral function: 
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lmi/j 5 (t) ~ 2f%M£6(t - M%) + "QCD continuum" 6 (i - t c 



(3) 



where the "QCD continuum comes from the discontinuity of the QCD diagrams, which is expected to give a 
good smearing of the different radial excitations []. The decay constant fjj is analogous to f„ = 92.32 MeV. 
However, in order to avoid some confusion, and for a more direct comparison with the lattice results, we 
shall abandon our favorite normalization, and adopt in the rest of the paper, the one: 



fo = V2f. 
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a normalization where = 130.56 MeV, and which the different experimental groups have also adopted; t c 
is the QCD continuum threshold, which is, like the sum rule variable t, an (a priori) arbitrary parameter. 
In this paper, we shall impose the t c and r stability criteria for extracting our optimal results The QCD 
expression of the correlator is well-known to two-loop accuracy (see e.g. [|| and the explicit expressions 
given in ||), in terms of the perturbative pole mass Mq, and including the non-perturbative condensates of 
dimensions less than or equal to six [J The sum rule reads: 



£(r) 




dt e~ tT 

M 2 87T 



H1 . xf i + !£)/(.) + £) 



R2s 



c A (o 4 



C 6 (0 6 )r e~ M Q T 



(5) 



where R2s is the new a^-term obtained semi-analytically in |0 and is available as a Mathematica package 
program Rvs.m. The other terms are: 



m = 



C 4 (0 4 ) 
C 6 (0 6 ) 



M 2 Q /t, 

^ + 2Li 2 (x) 

- log(l - x) 
-M Q (dd} + 



logxlog(l - x) - - log(l/x - 1) 



+ a;log(l/ir - 
(a s G 2 )/127r 



(x/(l - x))logx, 



1 



1 - -M^r g(da^G^d) 
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(6) 
(7) 



At the optimization scale, its effect is negligible, such that a more involved parametrization is not necessary. 
2 The corresponding t c value very roughly indicates the position of the next radial excitations. 

3 A different expression of the coefficient of the quark-gluon mixed condensate is given in jTj |, This change affects only 
slightly the result. We shall include the negligible contribution from the dimension six four-quark condensates. Notice that 
there is some discrepancy on the value of the four-quark coefficient in the literature. 
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The previous sum rules can be expressed in terms of the running mass toq(z^) through the perturbative 
three-loop relation p|, |l9|: 



M po i e = m{p 2 ) 



M 2 



Kq ' V 24 36 'V "' M- ' \ N 12 
, where, in the RHS, M po i e = M is the pole mass and: 

Kq = 17.1514 - 1.04137n + - V A f r = ^ 
w 3 ^ V Mr 

For < r < 1, A(r) can be approximated, within an accuracy of 1% by: 

A(r) ~ — r - 0.597r 2 + 0.230r 3 , 
8 

Throughout this paper we shall use the values of the parameters [pi Effl given in Table 1. 



(8) 



(9) 



(10) 



Table 1: Different sources of errors in the estimates of the decay constants (in MeV) and quark masses 
(in GeV). We have exagerately enlarged the error bars of different input in order to have conservative errors. 



Sources 



|A/ D |„- ic |A/ D | Mc |Am c | |AM C | |A/ B |^ \Mb\m„ |Am 6 | \AM b 



A 4 = (325 ± 43) MeV 
A 5 = (225 ± 30) MeV 
v = M C}b ± 1/2(M C)6 - m C)6 ) 
geom. estimate of a^-term 

t d = (1.2 ±0.2) GeV" 2 
t b = (0.35 ±0.05) GeV" 2 
6.0 < [GeV 2 ] < 9.5 
36. < tf [GeV 2 ] < 50 
(^^(l GeV)=-(229 ± 18) MeV 
(a s G 2 ) = (0.07 ±0.03) GeV 2 
Mq = (0.8 ±0.1) GeV 2 
a s (ipijj} 2 = (5.8 ±2.4) x lO" 4 GeV 6 
from our estimate of m c ^ or M c b 

Total 



7.4 


6.2 


0.03 


0.03 


















3.6 


3.0 


0.02 


0.02 


9.3 
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7.9 


8.2 






1.7 


2.0 






1.1 


1.1 


0.01 


0.01 
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5.7 


0.02 


0.02 


2.1 


2.8 


0.01 


0.01 


















2.4 


2.8 


0.05 


0.03 


8.8 


8.9 


0.01 


0.01 


5.2 


7.1 






1.8 


1.1 


0.005 


0.005 


0.9 


1. 






0.7 


0.7 


0.005 


0.01 


2.2 


2.1 


0.03 


0.04 


0.2 


0.3 






0.6 


1. 


0.01 


0.01 


10.9 


9.6 






15 


13 






20 


17 


0.04 


0.04 


23 


17 


0.06 


0.06 



We have used for the mixed condensate the parametrization: 

g(da^G^d) = M 2 (dd), (11) 
and we deduce the value of the QCD scale A from the value of a s (Mz) = (0.1184 ± 0.031) given in |2l|, q|. 

3 The D-meson channel 

We study in Fig. 1, the prediction of fn from the Laplace sum rules C and the one of Mo from ratio of 
moments 1Z for given value of the running charm quark mass m c (m c ). The influences of the choice of the 

4 It is clear that, for the non-perturbative terms which are known to leading order of perturbation theory, one can use either 
the running or the pole mass. However, we shall see that this distinction does not affect notably the present result. 
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continuum threshold t c and of the sum rule scale r are shown in details. Our optimal results correspond to 
the case where both stability in r and in t c are reached. However, for a more conservative estimate of the 
errors we allow deviations from the stability points, and we take: 

t c ~ (6 - 9.5) GeV 2 , r ~ (1.2 ± 0.2) GeV~ 2 , (12) 

and where the lowest value of t c corresponds to the beginning of the r-stability region. One can inspect that 
the dominant non-perturbative contribution is due to the dimension- four M c {dd) light quark condensate, 
and test that the OPE is not broken by high-dimension condensates at the optimization scale. However, the 
perturbative radiative corrections converge slowly, as the value of fn increases by 12% after the inclusion of 
the a s correction and the sum of the lowest order plus o^-correction increases by 21 % after the inclusion of 
the a 2 s term, indicating that the total amount of corrections of 21% is still a reasonnable correction despite 
the slow convergence of the perturbative series. However, as the radiative corrections are both positive, we 
expect that this slow convergence will not affect in a sensible way the final estimate. In order to improve the 
perturbative contributions, we have estimated the possible contribution of the af-term by assuming that its 
coefficient is the geometric sum of the a s and a 2 contributions. This effect is shown in Table 1, which is 
has a quite reasonnable value. A more precise answer on the higher order perturbative contribution needs 
an evaluation of the term which we hope to be available in the near future. In doing the analysis, one 
can also notice that the relative size of the perturbative corrections is smaller in the physical observable fn 
(perturbative+non-perturbative) than in the perturbative graph alone. This is due to fact that the role of 
the (iftip) condensate is important at the optimization scale, which then decreases the relative weight of the 
perturbative radiative corrections in the OPE. The estimate of Mn from the ratio of moments 1Z is less 
affected by radiative corrections, which tend to cancel each others due to the form of the sum rule. The 
behaviour of the optimized values of fn and Mn versus different values of m c (m c ) is given in Figs, la and 
lb. One can explicitly see in Figs, lc and Id that both Mn and fn are very sensitive to the change of 
rn c (m c ). This feature allows to have a good determination of the quark mass and then of fn, once the 
experimental value of Mn is used. Adding quadratically the different errors given in Table 1, we deduce the 
final estimate: 

m c (m c ) = (1.10 ± 0.04) GeV , f D = (201 ± 20) MeV , (13) 

where as mentioned previously, we have used the normalization f n — 130.56 MeV. These optimal values 
correspond to t c — 6.5 GeV 2 and t = 1.2 GcV~ 2 . The value of t c roughly corresponds to a radial excitation 
with a mass-splitting relative to the ground state mass of about M p which is phenomenologically acceptable. 
A similar analysis shown in Figs, le to lh is done for the pole mass. The discussions presented previously 
apply also here, including the one of the radiative corrections. We quote the final result: 

M c = (1.47 ± 0.04) GeV , f D = (208 ± 17) MeV , (14) 

where the error is slightly smaller here due to the absence of the subtraction scale uncertainties. For fn, we 
consider as a final estimate the mean value of the two predictions and taking the largest errors: 

fn = (204 ± 20) MeV . (15) 

Using our previous estimate of fnjfn = 1-15 ± 0.04 |lj|, we can also deduce: 

fn„ = (235 ± 24) MeV . (16) 

One can immediately compare the present predictions with the one obtained to order a s using the same 
procedure Q: 

fh c (m c ) = (1.08 ± 0.11) GeV , fn — (201 ± 15) MeV , (17) 

where one can notice a good agreement between the a s and a 2 results, though the error in [fl4f is smaller 
as the effect of the subtraction point v has not been taken into account. However, the agreement seems 
paradoxal in view of the fact that radiative corrections tend to increase the value of fn compared to the 
lower orders result. The different truncations of the expression of a s and of the relation between the pole 
and running mass also affect the absolute value of fn, which tend to compensate the increase due to the 
radiative corrections of the correlator. Therefore, a naive comparison becomes misleading. As one can see 
in Table 1, the main source of errors is due to the variations of the quark mass and to a lesser extent to the 
ones of v, A and (4>ip). The effect of t c on the result is relatively small from the value t c larger than 6 GeV 2 , 
where one starts to have a r stability. This result for the mass is also in agreement with the one from Mj/^ 
H §: 

m c (m c ) = (1.23tg:8£) GeV , (18) 
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but lower than the one from |p3| using non-relativistic Balmer formula for the cc bound state. One can 
cross-check that the two values of rfi c (m c ) and M c give the ratio: 

M c /fh c (m c ) ~ 1.33 , (19) 

which satisfies quite well the three-loop perturbative relation M c /fh c (m c ) = 1.33 obtained from the previous 
Eq. (||). This could be a non-trivial result if one has in mind that the quark pole mass definition can be 
affected by non-perturbative corrections not present in the standard SVZ-OPE. In particular, it may signal 
that 1/q 2 correction of the type discussed in || Q, if present, will only affect weakly the standard SVZ- 
phenomenology as observed explicitly in the light quark, gluonia and hybrid channels pif . Recent results to 
order a 2 [^2), using the analogous sum rule and using as input the pole mass value, gives fr> = (195 ± 20) 
MeV. The result is slightly lower than the result given here, though in agreement within the errors. We 
understand this slight difference as due to the lower value of the QCD continuum threshold used there, which 
corresponds to t c « 5.6 GeV 2 if one uses our pole mass value M c = 1.47 GeV. As shown in Fig. le, this t c 
value is on the boarder of the r stability region. It also indicates that the error introduced by the choice 
of t c could have been underestimated in this result. Fig. lg also indicates that fjj is quite sensitive in full 
QCD to the change of the M c value contrary to the remark given in Jl2| ]. Quenched and unquenched lattice 
results for fo are compiled in p5|. The quenched results range from (192 ± 18) MeV to (221 ± 17) MeV, in 



fair agreement with our results within the errors. The two available unquenched results 26 ^7j lead to the 
average: 

/jf = (220 ± 20) MeV , f§l = (254 ± 29) MeV , (20) 

which is again in agreement within the errors with our previous estimate. Finally, one can also compare the 
value of fr>^ with the experimental value [?]: 

/^ x p = (286 ± 60) MeV , (21) 

which agrees within la with our prediction. Improvements of our predictions for fjj s need an estimate of 
the ratio JdJ Id to order a 2 . 



4 The 5-meson channel 

We extend the previous analysis to the case of the B-meson, which again is an update of our previous work 
in |l4|, H |j| . The analysis is still similar to the one done in the D-channel, and is summarized in Fig 2 and in 
Table 1. Using the running &-quark mass, as a free parameter, we obtain at the optimization scale r = 0.375 
GeV" 2 and t c = 38 GeV 2 : 

rh b {m b ) = (4.05 ± 0.06) GeV , f B = (205 ± 23) MeV , (22) 

while using the pole mass as a free parameter, we get: 

AI h = (4.69 ± 0.06) GeV , f B = (200 ± 17) MeV , (23) 

from which we deduce the average: 

f B = (203 ± 23) MeV , (24) 

where we have taken the largest errors. One can again cross-check that the two values of rhbirrib) and Mb 
lead to 

M b /fhb(m b ) = 1.16 , (25) 

to be compared with 1.15 from the three-loop perturbative relation in Eq. (||), and might indirectly indicate 
the smallness of the l/<? 2 correction if any. Our result of fhb can be compared with our previous estimate 
from T-systems p2|| : 

fh b {m b ) = (4.23 ± 0.05) GeV , (26) 

and with similar values from recent estimates ^ |3^| . Our value of the perturbative pole mass is in agreement 
within the errors with the one in || ^2[ |3TJ] but is again lower than the one in p3| . Our value of f B is in fair 
agreement with the recent results f B = (206 ± 20) MeV obtained in |l2j using HQET sum rules and the one 
f B = (197 ± 23) MeV obtained in (l5) using the Laplace sum rule like in this work. The slight difference is 
due to the different appreciations of the continuum threshold t c and r stability regions in each papers. More 
specifically, errors related to the choice of t c at their choice of lower r-values appear to be underestimated 
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in these works. At such a choice of low r-values, the ground state contribution to the sum rule is smaller 
than in the present analysis. Unquenched lattice results [p7j p8| give the mean value: 

/]f = (198 ± 37) MeV , /s'/Zs* = 1.17 ± 0.03 , (27) 

where the largest error has been taken. These values agree with our previous determination in Eq. ( p4| ) and 
with our earlier estimate |l3[ : 

/ Bs // B = 1.16 ±0.05, (28) 

which has been confirmed from the recent analysis of jl5). Combining this SU(3) breaking ratio with our 
estimate of fs, one obtains: 

f Bs = (236 ± 30) MeV . (29) 
The extension of the previous analysis to the D* and B* channels is in progress. 

5 Summary 

We have updated our previous estimate of the quark masses and decay constants in |lj, |9[ || using the recent 
expression |l7j ] of a 2 a corrections for the heavy-light pseudoscalar correlators. Our results for the masses in 
Eqs. (jl|, [wf[22[ [H and for the decay constants in Eqs. |l(| ||, ||) confirm previous estimates to 



two-loops. The results for the running masses are: 

m c {m c ) = (1.10 ± 0.04) GeV , fh b {m b ) = (4.23 ± 0.05) GeV , Eqs @ ||) . 

The pole masses are: 

M c = (1.47 ± 0.04) GeV , M b = (4.69 ± 0.06) GeV , Eqs ([f| ||) . 

The decay constants are: 

f D = (204 ± 20) MeV , f B = (203 ± 23) MeV , Eqs (|| |§ . 

Using our SU(3) breaking prediction on fpjfp [Q, one also deduces: 

f Ds = (235 ± 24) MeV , f Bs = (236 ± 30) MeV , Eqs © ||) . 

The three-loop corrections tend to push the values of the decay constants to higher values restoring the 
slight discrepancy between the sum rules and recent unquenched lattice values. The resulting equality 
fu — fs confirm earlier findings from the sum rule |9) indicating large corrections to the 1/^/Mp heavy 
quark symmetry scaling law. Values of the quark masses obtained from the pseudoscalar sum rules are 
in agreement with the one from the quarkonia sum rules l| ^2|, [30|, but lower than the ones obtained in 
p3[ from non-relativistic Balmer formulae. The fitted values of the running and perturbative pole masses 
satisfy quite well their three-loop perturbative relation, which may indicate that 1 / g 2 -like terms & |24| have 
negligible effect in this channel like in the case of the light quark systems. 
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Figure 1: Laplace sum rule analysis of the decay constant /d, the running mass fh c (m c ) and the pole mass 
M c : a) fu versus the sum rule scale r at given m c (m c ) and for different values of the continuum threshold 
t c ; b) the same as a) but for Md- c) and d) effects of m c (m c ) on and fr>. The circle is the solution 
given by the data on M D ; e) to h) the same as a) to f) but for the pole mass. 
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Figure 2: The same as Fig. 1 but for the 6-quark and B-meson. 
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